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ABSTRACT 

We report the detection of 7-ray pulsations (> 0.1 GcV) from PSR J2229+6114 and 
PSR J1048— 5832, the latter having been detected as a low-significance pulsar by EGRET. Data 
in the 7-ray band were acquired by the Large Area Telescope aboard the Fermi Gamma-ray Space 
Telescope, while the radio rotational ephemerides used to fold the 7-ray light curves were obtained 
using the Green Bank Telescope, the Lovell telescope at Jodrell Bank, and the Parkes telescope. The 
two young radio pulsars, located within the error circles of the previously unidentified EGRET sources 
3EG J1048— 5840 and 3EG J2227+6122, present spin-down characteristics similar to the Vela pulsar. 
PSR J1048— 5832 shows two sharp peaks at phases 0.15 ± 0.01 and 0.57 ± 0.01 relative to the radio 
pulse confirming the EGRET light curve, while PSR J2229+6114 presents a very broad peak at phase 
0.49 ± 0.01. The 7-ray spectra above 0.1 GeV of both pulsars are fit with power laws having exponen- 
tial cutoffs near 3 GeV, leading to integral photon fluxes of (2.19 ± 0.22 ± 0.32) x 10~ 7 cm~ 2 s _1 for 
PSR J1048-5832 and (3.77±0.22±0.44) x Hr 7 cm- 2 s" 1 for PSR J2229+6114. The first uncertainty 
is statistical and the second is systematic. PSR J1048— 5832 is one of two LAT sources which were 
entangled together as 3EG J1048— 5840. These detections add to the growing number of young 7-ray 
pulsars that make up the dominant population of GeV 7-ray sources in the Galactic plane. 

Subject headings: gamma rays: observations — pulsars: general — pulsars: individual (PSR 
J1048— 5832, PSR J2229+6114) 
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1. INTRODUCTION 

The nature of unidentified high-energy 7-ray sources 
in the Galaxy was one of the major unanswered ques- 
tions at the end of the EGRET era. The third EGRET 
catalog contained 170 unidentified sources, 74 of which 
were at Galactic latitude \b\ < 10°(Hartman et al. 1999; 
Bhattacharya et al. 2003). Rotation-powered pulsars 
are believed to dominate the Galactic 7-ray source 



MS- 108, P. O. Box 1892, Houston, TX 77251, USA 

12 Istituto Nazionalc di Fisica Nuclcare, Sezione di Padova, 
1-35131 Padova, Italy 

13 Dipartimento di Fisica "G. Galilei" , Univcrsita di Padova, 
1-35131 Padova, Italy 

14 Department of Physics, Center for Cosmology and Astro- 
Particle Physics, The Ohio State University, Columbus, OH 43210, 
USA 

15 Istituto Nazionalc di Fisica Nucleare, Sezione di Perugia, 
1-06123 Perugia, Italy 

16 Dipartimento di Fisica, Univcrsita degli Studi di Perugia, 
1-06123 Perugia, Italy 

17 Dipartimento di Fisica "M. Merlin" dell'Universita e del 
Politccnico di Bari, 1-70126 Bari, Italy 

18 Istituto Nazionalc di Fisica Nucleare, Sezione di Bari, 70126 
Bari, Italy 

19 Laboratoire Leprincc-Ringuct, Ecole polytcchniquc, 
CNRS/IN2P3, Palaiseau, France 

20 Columbia Astrophysics Laboratory, Columbia University, 
New York, NY 10027, USA 

21 INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 
1-20133 Milano, Italy 

22 George Mason University, Fairfax, VA 22030, USA 

23 NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
USA 

24 Laboratoire de Physique et Chemie de l'Environnement, 
LPCE UMR 6115 CNRS, F-45071 Orleans Cedex 02, and 
Station de radioastronomic de Nangay, Obscrvatoire de Paris, 
CNRS/INSU, F-18330 Nancay, France 

25 Laboratoire de Physique Thcoriquc et Astroparticulcs, Uni- 
versity Montpellier 2, CNRS/IN2P3, Montpellier, France 

26 Department of Physics, Stockholm University, AlbaNova, 
SE-106 91 Stockholm, Sweden 

27 Royal Swedish Academy of Sciences Research Fellow, funded 
by a grant from the K. A. Wallenberg Foundation 

28 Dipartimento di Fisica, Universita di Udine and Istituto 
Nazionalc di Fisica Nucleare, Sezione di Trieste, Gruppo Collegato 
di Udine, 1-33100 Udine, Italy 

29 Universite de Bordeaux, Centre d'Etudes Nucleaires Bor- 
deaux Gradignan, UMR 5797, Gradignan, 33175, France 

30 CNRS/IN2P3, Centre d'Etudes Nucleaires Bordeaux Gradig- 
nan, UMR 5797, Gradignan, 33175, France 

31 Arecibo Observatory, Arecibo, Puerto Rico 00612, USA 

32 Department of Physical Sciences, Hiroshima University, 
Higashi-Hiroshima, Hiroshima 739-8526, Japan 

33 University of Maryland, College Park, MD 20742, USA 

34 University of Alabama in Huntsvillc, Huntsvillc, AL 35899, 
USA 

35 Australia Telescope National Facility, CSIRO, Epping NSW 
1710, Australia 

36 Department of Physics, Tokyo Institute of Technology, 
Meguro City, Tokyo 152-8551, Japan 

3 Waseda University, 1-104 Totsukamachi, Shinjuku-ku, Tokyo, 
169-8050, Japan 

38 Cosmic Radiation Laboratory, Institute of Physical and 
Chemical Research (RIKEN), Wako, Saitama 351-0198, Japan 

39 Department of Physics, University of Washington, Seattle, 
WA 98195-1560, USA 

40 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, BP 
44346, F-30128 Toulouse Cedex 4, France 

41 Jodrell Bank Centre for Astrophysics, School of Physics and 
Astronomy, The University of Manchester, M13 9PL, UK 

42 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigcl 
69, 53121 Bonn, Germany 

43 Center for Research and Exploration in Space Science and 
Technology (CRESST), NASA Goddard Space Flight Center, 
Greenbelt, MD 20771, USA 



population (e.g. Yadigaroglu & Romani (1995)), but 
their visibility is linked to their beam patterns. Soon 
after launch, the Fermi Gamma-ray Space Telescope 
began to unveil many 3EG sources, discovering the 
radio-quiet pulsar in the CTA1 supernova remnant 
associated with 3EG J0010+7309 (Abdo et al. 2008), 
detecting the radio-loud pulsar PSR J2021+3651 (asso- 
ciated with 3EG J2021+3716) (Abdo et al. 2009b), seen 
independently by AGILE (Halpcrn et al. 2008), as well 
as the radio pulsar PSR J1028— 5819, associated with 
3EG J1027-5817 (Abdo et al. 2009c) and new popula- 
tions of radio-quiet 7-ray pulsars, detectable using blind 
search techniques (Abdo et al. 2009f). For the pulsars 
detected in 7 rays, the bulk of the electromagnetic 
power output is in high energies. The 7-ray emission is 
thus crucial for understanding the emission mechanism 
which converts the rotational energy of the neutron star 
into electromagnetic radiation. The discovery of many 
new 7-ray pulsars will provide strong constraints on the 
location of the 7-ray emitting regions, whether above 
the polar caps (Daughcrty & Harding 1996), or far from 
the neutron star in the so-called "outer gaps" (Romani 
1996), or in the intermediate regions like the "slot gap" 
(Muslimov & Harding 2004) having "two-pole caustic" 
geometry (Dyks & Rudak 2003). In this paper we report 
the Fermi detection of the two pulsars PSR J1048— 5832 
and PSR J2229+6114, which have spin characteristics 
similar to other young pulsars typified by the Vela 
pulsar. Kramer et al. (2003) provide a discussion of 
"Vela-like" pulsars. 

PSR J1048-5832 (B1046-58) is located in the Carina 
region at low Galactic latitude (I = 287.42°, b = 0.58°). 
It was discovered during a 1.4 GHz Parkes survey of 
the Galactic plane and has a period P ^123.7 ms 
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(Johnston et al. 1992). It has a spin-down luminosity 
E = 4tt 2 I(P/P 3 ) of 2 x 10 36 ergs-\ for a moment 
of inertia / of 10 45 gcm 2 , a surface dipole magnetic 
field strength of 3.5 x 10 12 G, and a characteristic age 
t c = P/2P of 20kyr. High-resolution observations of 
the coincident ASCA source by the Chandra X-ray 
Observatory and XMM-Newton revealed an asymmetric 
pulsar wind nebula (PWN) of ~ 6" x 11", surrounding 
a point source coincident with the pulsar but so far no 
X-ray pulsations were detected (Gonzales et al. 2006). 
PSR J1048— 5832 has been proposed as a counterpart 
of the steady EGRET source 3EG J1048-5840 (Fierro 
1995; Pivovaroff et al. 2000; Nolan et al. 2003), as 
suggested by positional coincidence, spectral, and 
energetic properties. Detailed analysis of the EGRET 
data found possible 7-ray pulsations at E >400 MeV, 
a double-peaked light curve with ~ 0.4 peak phase 
separation (Fig.l, Kaspi et al. 2000). An HI distance 
determination for the pulsar yields between 2.5 and 6.6 
kpc (Johnston et al. 1996). The NE2001 model (Cordes 
& Lazio 2002) assigns a distance of 2.7 kpc based partly 
on the HI distance determination. For this paper we 
adopt 3 kpc as the distance to the pulsar. 

PSR J2229+6114 is located at (l,b) = (106.6°,2.9°) 
within the error box of the EGRET source 
3EG J2227+6122 (Hartman et al. 1999). Detected 
as a compact X-ray source by ROSAT and ASCA 
observations of the EGRET error box, it was later 
discovered to be a radio and X-ray pulsar with a period 
of P = 51.6 ms (Halpcrn et al. 2001b). The radio pulse 
profile shows a single sharp peak, while the X-ray light 
curve at 0.8 - 10 keV consists of two peaks, separated 
by A(f> = 0.5. AGILE recently reported the discovery of 
7-ray pulsations above 100 MeV (Pellizzoni et al. 2009). 
The pulsar is as young as the Vela pulsar (characteristic 
age r c = lOkyr), as energetic (E = 2.2 x 10 37 crgs _1 ), 
and is evidently the energy source of the "Boomerang" 
arc-shaped PWN G106.65+2.96, suggested to be part of 
the supernova remnant (SNR) G106.3+2.7 discovered 
by Joncas & Higgs (1990). Recently, the PWN has been 
detected at TeV energies by MILAGRO (Abdo et al. 
2009g). Studies of the radial velocities of both neutral 
hydrogen and molecular material place the system at 
~ 800 pc (Kothes et al. 2001), while Halpern et al. 
(2001a) suggest a distance of 3 kpc estimated from its 
X-ray absorption. The pulsar DM, used in conjunction 
with the NE2001 model, yields a distance of 7.5 kpc, 
significantly above all other estimates. For this paper 
we again adopt a distance of 3 kpc. 

2. OBSERVATIONS 

2.1. Gamma-Ray Observations 

The Large Area Telescope (LAT) aboard Fermi is 
an electron-positron pair conversion telescope sensitive 
to 7-rays of energies 0.02 to >300 GeV. The LAT is 
made of a high-resolution silicon microstrip tracker, a Csl 
hodoscopic electromagnetic calorimeter, and segmented 
plastic scintillators to reject the background of charged 
particles (Atwood et al. 2009). Compared with its pre- 
decessor EGRET, the LAT has a larger effective area 
(—8000 cm 2 on-axis for E > 1 GeV), improved angu- 
lar resolution (# 6 8 ~ 0.5° at 1 GeV for events in the 



front section of the tracker), and a higher sensitivity 
(~3x lO^cm^s- 1 ) 60 . The large field of view (-2.4 
sr) allows the LAT to observe the full sky in survey mode 
every 3 hours. The LAT timing is derived from a GPS 
clock on the spacecraft and 7-rays are hardware time- 
stamped to an accuracy significantly better than 1 /is 
(Abdo et al. 2009h). The LAT software tools for pulsars 
have been shown to be accurate to a few fis (Smith et al. 
2008). 

In this paper, the data used for the spectral analysis 
were collected during Fermi's first-year all-sky survey, 
beginning 2008 August 4. For the timing analysis, we 
added data collected during the commissioning phase ob- 
servations from 2008 June 30 to 2008 August 3, that in- 
cluded pointed observations of the Vela pulsar and other 
targets. The data for PSR J1048-5832 end 2009 April 
10, while those for PSR J2229+6114 end 2009 March 23. 
Only 7 rays in the "Diffuse" class events (the tightest 
background rejection) were selected. In addition, we ex- 
cluded those coming from zenith angles > 105°, where 7- 
rays resulting from cosmic-ray interactions in the Earth's 
atmosphere produce an excessive background contamina- 
tion. 

2.2. Radio Observations 

2.2.1. PSR J '1048- 5832 

PSR J1048-5832 is observed in the radio band at ap- 
proximately monthly intervals since early 2007 at the 64- 
m Parkes radio telescope in Australia (Manchester 2008). 
A typical observation is of 2 min duration at a frequency 
of 1.4 GHz with occasional additional observations at 0.7 
and 3.1 GHz. Full details of the observing and data anal- 
ysis can be found in Wcltevrede et al. (2009a). The 
pulsar is known to have glitched in the past (Wang et 
al. 2000) and indeed glitched just prior to the launch of 

Fermi (Weltevrede et al. 2009a). Like many high E pul- 
sars, PSR J1048— 5832 is highly polarized in the radio 
band over a wide frequency range (Karastergiou et al. 
2005; Johnston et al. 2006). 

The timing solution uses 20 pulse times of arrival 
(TOAs) and was derived using the pulsar timing software 
TEMP02 (Hobbs et al. 2006). The timing solution fits 
for the pulsar's spin frequency and frequency derivative 
and also whitens the data using the so-called fitwaves al- 
gorithm within TEMP02. The resulting rms of 0.287 ms 
is a substantial improvement on the prc-whitencd solu- 
tion. Note that the extra fitwaves parameters are not 
supported by the standard Fermi software tools but were 
included in this analysis. Measurements of the dispersion 
delay across the 1369 MHz band with a bandwidth of 256 
MHz lead to a dispersion measure (DM) of 128.822 ± 
0.008 pccm~ 3 , with no indication that the DM is vary- 
ing over time. This DM is used to correct the radio TOAs 
to infinite frequency for phasing the 7-ray photon TOAs. 

2.2.2. PSR J2229+6U4 

PSR J2229+6114 is being observed at the NRAO 
Green Bank Telescope (GBT) (Kaplan et al. 2005) and 

60 This value refers to a steady source after one-year sky survey, 
assuming a high-latitude diffuse flux of 1.5 X 10~ 5 cm~ 2 s _1 sr -1 
( >0.1 GeV and a photon spectral index of -2.1 with no spectral 
cutoff). 
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the Lovell telescope at Jodrcll Bank (Hobbs et al. 2004). 
The rotational ephemeris used here to fold 7-ray pho- 
tons is based on TOAs obtained from both telescopes 
between 2008 June 17 and 2009 March 23. There are 
25 such TOAs from GBT with average uncertainty of 
0.2 ms, each from a 5 minute observation mainly at 
a central frequency of 2.0 GHz. The 44 Jodrell Bank 
TOAs have an average uncertainty of 0.3 ms, obtained 
from individual 30 minute observations at 1.4 GHz. In 
our timing fits with TEMPO 61 , we fixed the position 
at that known from X-ray observations (Halpern et al. 
2001b). PSR J2229+6114 shows some timing noise over 
the 9 month interval, requiring a fit to rotation frequency 
v = l/P and its first two derivatives. In addition, on 
MJD 54782.6 ± 0.5 a small glitch occurred, with frac- 
tional frequency step of Avjv = (4.08 ±0.06) x 10~ 9 and 
Av/i> = (2.0 ± 0.4) x 10~ 4 . The post-fit rms is 0.24ms. 
We have measured the DM with separate sets of GBT ob- 
servations at three widely spaced frequencies, obtaining 
DM = (204.97 ± 0.02) pc cm" 3 . This is used to correct 
2 GHz arrival times to infinite frequency, for comparison 
with the 7-ray profile. 

3. ANALYSIS AND RESULTS 

The events were analyzed using the standard software 
package Science Tools 62 (ST) for the Fermi LAT data 
analysis and TEMP02. The timing parameters used in 
this work will be made available on the servers of the 
Fermi Science Support Center 63 . 

3.1. Pulse profiles of PSR J1048-5832 

For the timing analysis of PSR J1048-5832, we se- 
lected 7-rays with energy >0.1 GeV within a radius of 
1?0 around the radio pulsar position. Then, we ap- 
plied an energy-dependent angular radius cut O^s ^ 
0.8 x Eq°J 5 de grees, keeping all the photons included in 
a radius of 0?35. This selection approximates the LAT 
point spread function (PSF) and maximizes the signal- 
to-noise ratio over a broad energy range. We corrected 
photon arrival times to the Solar System barycenter us- 
ing the JPL DE405 Solar System ephemeris (Standish 
1998). The event times were then folded with the radio 
period using the Parkes ephemeris. The bin-independent 
-ff-test (De Jager et al. 1989) results in a probability 
< 4 x 10 -8 that the modulation would have occurred by 
chance. This value is more than 4 orders of magnitude 
smaller than the previous EGRET results (Kaspi et al. 
2000), establishing firmly this source as a 7-ray pulsar. 

Figure 1 (top panel) shows the resulting 50-bin 7-ray 
phase histogram above 0.1 GeV. The two peaks appear 
asymmetric, with a slow rise and a fast fall. We fit 
each peak with two half-Lorentzian functions, i.e. with 
different widths for the leading and trailing sides. PI 
((f) = 0.05—0.17) is sharper with a full width at half maxi- 
mum (FWHM) of 0.06 ±0.01, while P2 (0 = 0.45-0.65) 
is a bit broader with an FWHM of 0.10 ± 0.02. The 
1.4 GHz radio profile is shown at the bottom panel for 
comparison. The 7-ray light curve consists of two peaks, 
with PI at phase 0.15 ± 0.01 ± 0.0001 and P2 at phase 
0.57± 0.01 ± 0.0001, leading to a phase separation A(f> of 

61 http://www.atnf.csiro.au/research/pulsar/tempo. 

62 http://fcrmi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview. 

63 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/cphems/ 



0.42±0.01±0.0001. The errors are respectively the 7-ray 
fit uncertainty and that caused by the DM uncertainty. 
These results are in agreement with those found by Kaspi 
et al. (2000). We observe structure related to a shoul- 
der or a "bridge" region between 0.17 and 0.30 in phase, 
which trails the first 7-ray peak. Note that the profile 
is very similar to the Vela light curve, which consists of 
both two peaks separated by 0.43 in phase and a bridge 
region (Abdo et al. 2009a). We also defined the "off- 
pulse region" as the pulse minimum (<p — 0.7— 1.05). To 
check this assumption, we estimated the background rep- 
resented by the dashed line (73 counts/bin) from a ring 
with 1°< 9 <2° surrounding the source. Nearby sources 
are removed, and we normalized to the same phase space 
as our selection. The result is in good agreement with 
the off-pulse region. As a consequence, the total number 
of pulse photons from the pulsar is estimated at 933 ±93, 
with a background contribution of 3654 ± 60 events. 

Figure 1 (middle panels) shows the 50-bin 7-ray phase 
histograms of PSR J1048— 5832 in four energy bands (0.1 
- 0.3 GeV, 0.3 - lGeV, 1 - 3 GeV, >3GeV). Evolution 
in the light curve shape with energy is visible, although 
more data are needed to constrain the peak widths. No- 
tably, below 0.3 GeV the first peak seems wider than at 
high energies. This feature could be explained by the 
contamination of a LAT source (0FGL 1045.6—5937) less 
than 1° from the pulsar. We also looked for an evolution 
of the ratio P1/P2 with energy, as seen by EGRET for 
Vela, Crab, Geminga, PSR B1951±32 (Thompson et al. 
2004) and seen by the Fermi LAT for Vela (Abdo et al. 
2009a) and PSR J0205+6449 (Abdo et al. 2009e). In 
each energy band, we calculated the peak height with 
respect to the background level as estimated from the off 
pulse interval. The ratio shows no variation with a confi- 
dence level of ~69%. This is not a particularly stringent 
bound, and we expect to accumulate much more data 
to eventually detect a possible variation of P1/P2 with 
energy. Finally, note that between 1 and 3 GeV, the 
trailing shoulder PI ((f) — 0.17 - 0.33) is significant and 
there is still evidence of it for E > 3 GeV, whereas above 
3 GeV the two peaks are still observed with the highest 
energy photon detected in PI at 19 GeV. 

3.2. Pulse profiles of PSR J2229+6114 

For the timing analysis of PSR J2229+6114, we used 
the same energy- dependent selection criteria as for PSR 
J 1048- 5832. 7-ray TOAs were folded with the radio 
period using the GBT and Lovell telescope ephemeris. 

Figure 2 (top panel) shows the resulting 50-bin his- 
togram of folded counts above 0.1 GeV compared with 
the phase-aligned radio pulse profile (bottom panel). The 
profile shows a single, asymmetric peak ((f) = 0.15 — 0.65), 
that has been fit by a two half-Lorentzian function. The 
fit places the peak at 0.49 ± 0.01 ± 0.001 with FWHM of 
0.23 ±0.03. We estimated the background level from a 1 

2° ring around the pulsar during the 0.65 - 0.15 pulse 
minimum. This is represented by the dashed line (48 
counts/bin). The result is consistent with the off-pulse 
region. As a consequence, the total of pulsed photons 
from the pulsar is estimated at 1365 ± 79, with a back- 
ground contribution of 2431 ± 49 events. 

To examine the energy-dependent trend of the 7- 
, m ray pulse profile (Figure 2, middle panels), phase his- 
tograms are plotted over 4 energy intervals (0.1—0.3 GeV, 
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Fig. 1. — Top panel: Light curve of PSR J1048— 5832 above 0.1 GeV, shown over two pulse periods with 50 bins per period (P = 123.7 ms). 
The dashed line shows the background level, as estimated from an annulus surrounding the pulsar position during the off-pulse phase (73 
counts/bin). Four following panels: Energy-dependent phase histograms for PSR J1048— 5832 in the four indicated energy ranges, each 
displayed with 50 bins per pulse period. Second panel from bottom: Folded EGRET light curve for energies above 400 MeV (Kaspi ct al. 
2000). Bottom panel: Radio pulse profile from the Parkes Telescope at a center frequency of 1.4 GHz with 1024 phase bins. 
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Fig. 2. — Top panel: Light curve of PSR J2229+6114 above 0.1 GeV, shown over two pulse periods with 50 bins per period (P = 51.6 ms). 
The dashed line shows the background level, as estimated from an annulus surrounding the pulsar position during the off-pulse phase (48 
counts/bin). Four following panels: Energy dependent phase histograms for PSR J2229+6114 in the four indicated energy ranges, each 
displayed with 50 bins per pulse period. Second panel from bottom: Light curve in the 1—10 keV band from the XMM pn CCD in small 
window mode. The instrumental time resolution is 5.7ms (2.2 phase bins). Phase alignment with respect to the radio pulse, as described 
in the text, is accurate to ?s 0.4 phase bins of this 20-bin light curve. Bottom panel: Radio pulse profile from Green Bank Telescope at a 
center frequency of 2 GHz with 128 phase bins. 
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0.3 — 1 GeV, 1 — 3 GeV and >3 GeV), showing a possible 
drift of the peak. Between 0.1 - 0.3 GeV, the peak is 
offset from the radio pulse by 0.51 ± 0.02 according to a 
two half-Lorentzian fit, while the offsets for 0.3 - 1 GeV, 
1 - 3 GeV and > 3 GeV are 0.48 ± 0.01, 0.49 ± 0.01 and 
0.45 ±0.01, respectively. The peak positions are compat- 
ible between 0.1-3 GeV, and a slight mis- alignment ap- 
pears above 3 GeV. However, the data cannot constrain 
strongly the pulse phase dependence, even if the mis- 
alignment between the X-ray and 7-ray profile suggests 
a spectral energy dependence of the high energy light 
curve. Finally, we note that the highest energy photon 
is 10.8 GeV at phase ~ 0.30. 

In Figure 2 we also show the 1-10 keV X-ray profile of 
PSR J2229+6114 obtained from an XMM observation on 
2002 June 15 (MJD 52440) with effective exposure time 
of 20 ks. The data were folded using a contemporaneous 
ephemeris based on Jodrell Bank and GBT observations. 
The highest peak of the X-ray profile lags the radio pulse 
by 4> — 0.17 ± 0.02. Analysis of an RXTE observation 
from MJD 52250 yields a consistent radio-X-ray offset. 
There is no energy dependence of the X-ray pulse shape 
within the 1-10 keV range, while the sharpness of the 
peaks as well as their spectral shape indicates that the 
emission is predominantly non-thermal. 

3.3. Phase-averaged spectra and flux 

In order to obtain the phase-averaged spectrum for 
PSR J1048-5832 and PSR J2229+6114, a maximum 
likelihood spectral analysis (Mattox et al. 1996) was per- 
formed, using the LAT tool 'gtlike' 3 . We used the "Pass 
6 v3" instrument response functions (IRFs) , which are a 
post-launch update to address 7-ray detection inefficien- 
cies that are correlated with trigger rate. The systematic 
errors on the effective area are < 5% near 1 GeV, 10% 
below 0.1 GeV and 20% over 10 GeV. 

3.3.1. Spectrum of PSR Jl 048- 5832 

A circular region around 15° from the source was mod- 
eled including nearby LAT sources around the source po- 
sition. The Galactic diffuse background was accounted 
for by using maps based on the GALPROP model called 
54_59Xvarh7S (Strong et al. 2004a,b). We modeled the 
isotropic background through a tabulated spectrum de- 
rived from a fit to LAT data at high galactic latitude, 
freezing the Galactic model and including the detected 
point sources within a radius of 15°. We fit the spectrum 
of PSR J1048— 5832 with a power law with an exponen- 
tial cutoff between 0.1 and 0.7 in phase, that can be 
described by the equation: 

% = N E- T e- E ' E ' cm-V^cV- 1 (1) 
dE 

with E in GeV, the term N Q = (5.9 ± 0.3 ± 0.1) x 
10~ 8 cm -2 s -1 GeV -1 , a spectral index Y = 1.38±0.06± 
0.12 and a cut-off energy E c = 2.'3±° Q l ± 0.3 GeV. The 
first uncertainty is statistical, while the second is system- 
atic. From this result, we obtained an integral photon 
flux in the range 0.1-100 GeV of (2.19 ± 0.22 ± 0.32) x 
10 _7 cm _2 s _1 , which is about one third of the flux of 
the EGRET source 3EG J 1048- 5840. Figure 3 shows 
both the overall fit between 0.1 and 20 GeV (solid line), 



and the spectral points obtained for 6 logarithmically- 
spaced energy bins and performing spectral analysis in 
each interval, assuming a power law shape for the source. 
To check the assumption of a cut-off energy in the spec- 
trum, we also modeled the pulsar with a simple power 
law of the form dF/dE = N a E~ r . The chance probabil- 
ity to incorrectly reject the hypothesis of a pure power 
law spectrum is ~ 10 a. 

The pulsar is listed in the Bright Source List of the 
Fermi LAT (Abdo et al. 2009d) as 0FGL J1047.6-5834, 
which is located at (RA,Dec) = (161.922°, -58.577°) 
with a 95% confidence level radius of 0.138°. There 
is another LAT point source located ~1° away, 
0FGL J1045.6-5937. We modeled this unidentified 
source with a simple power law. The best fit result gives 
a spectral index 2.2 ± 0.1 ± 0.1 and an integral photon 
flux of (4.49 ± 0.40 ± 0.80) x 10~ 7 cm" 2 s" 1 . The sum 
of the fluxes of this source plus PSR J 1048- 5832 is ~ 
6.7 xl0~ 7 cm _2 s _1 which is consistent with the flux of 
3EG J1048-5840 of (6.2 ± 0.7) x rO^cni^s" 1 . The 
EGRET source 3EG J1048-5840 as well as the COS-B 
source 2CG 288-00 (Swanenburg et al. 1981) were ap- 
parently made up of these two sources, which have now 
been resolved by the Fermi LAT. 

3.3.2. Spectrum of PSR J2229+6449 

Initially, a 15 degree circular region around the 
source position and the same Galactic diffuse back- 
ground used for PSR J1048-5832 were modeled to fit 
PSR J2229+6114. However, some structures around the 
object were not taken into account in the Galactic model, 
and overestimate the flux of the pulse regions. We finally 
adopted a model for the Galactic diffuse emission based 
on six galactocentric 'ring' maps of N(H I) and W(CO) 
and on the spatial distribution of the inverse Compton 
intensity modeled by GALPROP. This intensity as well 
as the gamma-ray emissivities per ring were adjusted to 
maximize agreement with the observations, taking into 
account detected point sources of gamma rays. This ap- 
proach is similar to that described by Casandjian & Gre- 
nier (2008) for modeling EGRET data, and is being used 
by the LAT team to fit the model of Galactic diffuse emis- 
sion for the first public release. 

PSR J2229+6114, referenced as the LAT point source 
0FGL J2229.0+6114 in the Bright Source List, was mod- 
eled by a power law with a simple exponential cut- 
off (Eq. 1), between 0.15 and 0.65 in phase. Fig- 
ure 4 (solid line) shows the phase-averaged spectral en- 
ergy distribution from the likelihood fit, with Nq — 
(5.2±0.4±0.1)x 10- 8 cm- 2 s- 1 GcV- 1 , a spectral index 
r = 1.85 ±0.06 ±0.10 and a cut-off energy E c = 3.6^ ± 
0.6 GeV. From this, we estimated an integral photon flux 
of (3.77 ± 0.22 ± 0.44) x lO^cm" 2 s" 1 . This value is 
10% lower than the flux of (4.13 ±0.61) x 10~ 7 cnr 2 s" 1 
of the EGRET source 3EG J2227±6122 obtained for a 
power law with a spectral index of 2.24 ± 0.14 and larger 
than the flux of (2.6 ± 0.4) x 10~ 7 cm _2 s _1 measured 
by AGILE (Pcllizzoni et al. 2009). We also fit the 
LAT data to a simple power-law model, yielding a flux 
of (4.54 ± 0.14) x 10~ 7 cm~ 2 s~ 1 and a spectral index 
of 2.25 ± 0.02. Note that the spectral model using an 
exponential cut-off is better constrained with a differ- 
ence between the log likelihoods of ~ 9a, rejecting the 
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Fig. 3. — Spectral energy distribution for PSR J1048-5832 as fit by "gtlike" assuming a power-law with an exponential cutoff (solid line). 
The spectral points were obtained for 6 logarithmically-spaced energy bins and performing spectral analysis in each interval, assuming a 
power law shape for the source. The errors bars are statistical only. 



power-law hypothesis. For both PSR J1048-5832 and 
PSR J2229+6114 the statistics available is not enough 
to significantly rule out a super-exponential cutoff in fa- 
vor of a simple exponential cutoff. 

As a first search for unpulscd emission from the 
Boomerang nebula, we fitted a point-source to the off- 
pulse data at the radio pulsar position in the energy 
band 0.2 - 100 GeV. No signal was observed from the 
PWN. After scaling to the full pulse phase, we de- 
rived a 95% confidence level upper limit on the flux of 
4.0 x 10- 8 cm _2 s _1 . 

4. DfSCUSSfON 

The 7-ray light curves of PSR J1048-5832 and 
PSR J2229+6114 cover a wide range in phase, suggest- 
ing that the 7-ray beams cover a large solid angle. This 
in turn seems to favor the outer magnetospheric emis- 
sion models, in particular the outer gap (OG, Romani 
1996) and the slot gap (SG, Muslimov & Harding 2004) 
models. 

The 7-ray luminosity is given by 

£7 =^fUaXE)F E ,obsd 2 1 (2) 

where FE,obs is the observed energy flux at the Earth line 
of sight (at angle Qe to the rotation axis), d is the pulsar 
distance, and fa(a, Ce) is the beaming correction factor 
that depends on the geometry of the emission pattern. 
The factor fa(a, Ce) 1S a function of the pulsar magnetic 
inclination a, and is model sensitive. It is given by (Wat- 
ters et al. 2009): 



/ F 7 (a;C» sin(()d(d<f> 
2fF 1 (a;( E , 



(3) 



where Fj(a; C, 4>) is the radiated flux as a function of the 
viewing angle C an d the pulsar phase 4>. In the ratio 



fa, the numerator is the total emission over the full sky, 
and the denominator is the expected phase-averaged flux 
for the light curve seen from Earth. For the polar cap 
model (Daugherty & Harding 1996), the 7-ray emission 
originates at a few stellar radii from the surface, implying 
an emission with a small solid angle, that is, fa <C 1. For 
both OG and SG models, where the emission is far away 
from the neutron star, the resulting fa values can be near 
unity, and even exceed 1 when the dominant sampling of 
the flux F 7 arises for viewing angles £ quite disparate 
from 

For PSR J1048-5832, the observed integral energy flux 
obtained by integrating Equation 1 times the energy is 
(19. 4± 1.0±3.1) x 10~ n ergcm~ 2 s _1 , leading to a 7-ray 
luminosity of 2.1 x 10 35 fa (d/3 kpc) 2 and an efficiency 
r) = L-t/E = 0.10 fa (d/3 kpc) 2 . 

The peak separation in the 7-ray light curve and 
the 7-ray efficiency are useful to place constraints in 
the now- favored outer magnetospheric emission models. 
Using the 7-ray light curve "Atlas" of Watters et al. 
(2009) we can estimate an allowed range a ~ 60°- 85° , 
C ~ 70°- 80° , and fa ~ 0.7 - 1.1 for the OG model, 
and a ~ 50°- 75° , C ~ 50°- 75° , and fa ~ 1.1 for 
the SG (Two-Pole Caustic) model. The variation of the 
radio polarization position angle for a pulsar constrains 
the impact angle (3e = en - Ce- Fitting the rotating 
vector model (Radhakrishnan & Cooke 1969) to the 
PSR J1048— 5832 radio data, a value of /3g smaller 
than 10° has been derived (Karastergiou et al. 2005; 
Weltevrede & Johnston 2008a) . Both theoretical models 
thus have a good range of possible solutions, but these 
estimates assume the efficiency relation rj ~ (1Q 33 / E) - 5 
of Watters et al. (2009), which gives r\ ~ 0.02, about a 
factor of five smaller than that derived for this pulsar 



9 



T > > > I | 




_l I I I I I I I I I I I I I I I I I I 

1 0" 1 1 10 

Energy (GeV) 

Fig. 4. — Spectral energy distribution for PSR J2229+6114 as fit by "gtlike" assuming a power-law with an exponential cutoff (solid line). 
The spectral points were obtained for 6 logarithmically-spaced energy bins and performing spectral analysis in each interval, assuming a 
power law shape for the source. The errors bars are statistical only. 



from its inferred luminosity. Assuming our measured 
quantity 77 = 0.10 the allowed range remains about the 
same for the Two-Pole Caustic (TPC) model and the 
range of a shifts to ~ 70° - 90° for the OG model. The 
phase lag between the radio pulse and the first 7-ray 
pulse is in agreement with both OG and TPC or SG 
models. 

For PSR J2229+6114, the observed integrated energy 
flux is (23.7 ± 0.7 ± 2.5) x lO" 11 ergcuT 2 s _1 . Thus, we 
estimate a 7-ray luminosity L 1 = 2.6x 10 35 /n(rf/3 kpc) 2 
erg s _1 , and deduce an efficiency ry 7 = 0.011 fa 
(d/3 kpc) 2 , which is a factor 1.6 larger than the effi- 
ciency estimated by the relation 77 ~ (10 33 /-E) - 5 . Note 
that with the estimated SNR distance of 0.8 kpc the ef- 
ficiency decreases to 0.001 fa, while with the distance 
derived from the DM (7.5 kpc) the efficiency increases 
to 0.07 fa, emphasizing the importance of the distance 
determination. 

As for PSR J1048— 5832, we can compare light curves 
with the geometrical models of Watters et al. (2009) 
to derive constraints on the geometry. For outer mag- 
netosphere models with ?7 7 ~ 0.01, both TPC and OG 
geometries can deliver single 7-ray pulses at a range of 
angles. Single pulse solutions appear for a ~ 20° — 55°, 
C - 25° - 50° for the TPC/SG geometry and over the 
range a ~ 45° - 80° , C ~ 35° - 70° and fa ~ 0.47 - 1.08 
for the OG model. For PSR J2229+6114, we have an 
additional geometrical constraint from modeling of the 
X-ray PWN torus surrounding the pulsar, as measured 
by Chandra; Ng & Romani (2004, 2008) fit these data 
to determine a viewing angle £ = 46° ± 2° ± 6° . This is 
in the overlap range consistent with both the TPC and 
the OG single pulse solutions. Finally, although there 



is some uncertainty in determining the precise phase 
of the magnetic axis from the radio pulse, the lag of 
the strong 7-ray peak from the radio peak (0 w 0.50) 
provides additional information. Examining the sample 
light curves in Figures 9 & 10 of Watters et al. (2009) we 
find that single peaks at this phase are found only when 
they can be identified as the normal 'P2' component, 
with a 'bridge' of emission to earlier phases and 'PI' 
faint or missing. Such light curves appear only for small 
efficiency 7/ < 0.03 and only for a modest range of angles: 
a ~ 50° - 70°, C ~ 40° - 50° (TPC/SG, with PI faint) 
or a > 45°, C = 45° - 50° (OG, with PI absent). This is 
precisely the £ range required by the X-ray torus fitting; 
including the small (3e constraint from the detection of 
radio emission, we see that the radio, X-ray and 7-ray 
data are all consistent with a = 55° ± 5°, Q — 45 ± 5°, 
implying fa 1 for both models. Thus low 7-ray 
efficiency, high altitude emission and a well constrained 
viewing geometry together provide a consistent picture 
for gamma-ray light curve shape. On the other hand, 
the gamma-ray light curve does not seem to match the 
the clearly double-peaked X-ray light curve. However, 
we have argued that the gamma-ray light curve is 
intrinsically double-peaked, with the first peak missing 
or very weak, and the appearance of a possible first peak 
in the > 100 MeV an d 0.1 - 0.3 gamma-ray light curves 
(signal/-y/background = 2.5 a) that is around the phase 
of the first X-ray peak is intriguing. But the second 
gamma-ray peak is not in phase with the second X-ray 
peak. As the radiation mechanisms contributing to the 
emission in X-ray and LAT wavebands may differ, and 
at a given magnetospheric colatitude, the anisotropy 
of each emission mechanism is both energy- dependent 
and altitude-dependent, one does not expect the pulse 
profiles to be achromatic. Therefore, more detailed 
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physical radiation models are required to understand 
the energy dependence of the LAT light curve and to 
explain the lower energy X-ray light curve. 

This evidence against low altitude emission in these 
pulsars can also be supplemented by constraints of a sep- 
arate physical origin. The observation of photons out to 
beyond 5 GcV precludes any dominant action of mag- 
netic pair creation 7 — ► e + e~ in the emission region 
at energies below this. Accordingly, the maximum ob- 
served energy, e max GeV, provides a lower bound to the 
altitudes of emission, since it must lie below any thresh- 
old energy for a super-exponential 7-B pair production 
turnover. Such constraints have been obtained for Fermi 
detections of the Vela pulsar (Abdo et al. 2009a), and 
PSR J1028-5819 (Abdo et al. 2009e), indicating that the 
super-GeV emission must originate above at least 2.2 and 
2.1 stellar radii, respectively. Using a standard polar 
cap model estimate for the minimum emission height of 
r > {e max B 12 /1.76GeV) 2 / 7 P' 1 / 7 R* (e.g. inverting Eq. 
[1] of Baring (2004)), for a surface polar field strength 
of 10 12 Bi 2 G, the PSR J1048-5832 spin-down parame- 
ters (P — 0.1237s, B12 — 3.5) together with the max- 
imum observed energy e max ~ 9 GeV in Figure 3 yield 
r > 3.LR*. This drops to r > 2.1R* if E c ~ 2.3 GeV 
is deployed when defining this bound (as was the case 
for the aforementioned Fermi pulsar detections). Simi- 
larly, for PSR J2229+6114 with P = 0.0516 s, B 12 = 6.4 
and the maximum observed energy e max ~ 8 GeV in Fig- 
ure 4, the altitude constraint is r > 4.0i?». For either 
case, clearly these bounds preclude emission very near 
the stellar surface, adding to the advocacy for a slot gap 
or outer gap acceleration locale for the emission in both 
pulsars. 

5. CONCLUSION 

Although PSR J1048-5832 and PSR J2229+6114 arc 
both Vela-like in age and spin characteristics, their light 
curves and derived emission geometries are quite differ- 
ent. Table 1 summarises the main quantities measured 
for both pulsars. The double-peaked light curve of PSR 
J 1048— 5832 is nearly identical to that of Vela, whereas 
its derived efficiency is a factor of 10 larger than that of 
Vela (Abdo et al. 2009a), adopting /n = 1 and d — 3kpc. 
On the contrary, the 7-ray efficiency of J2229+6114 is 
very similar to that of the Vela pulsar, but the pulsar 
shows a single, large peak similar to PSR J0357+32 dis- 
covered by searching for pulsations at the positions of 
bright 7-ray sources (Abdo et al. 2009d). Note that the 
efficiency of PSR J2229+6114 would be smaller at the 
distance of about lkpc that some authors suggest. 

With the growing number of detected 7-ray pulsars, we 
are beginning to sample a wider variety of emission and 
viewing geometries, and pulsar ages. The range of light 
curve morphologies should allow improved constraints on 
high-energy emission models and a better understanding 



of the pulsar magnetospheric structure and acceleration 
process. For example, while many young pulsars, like 
J1048— 5832, show Vela-type light curves, a small num- 
ber are similar to J2229+6114, with a strong P2 compo- 
nent and a weak or absent PI (Weltevrede et al. 2009b). 
Mapping the angle range over which PI is missing, espe- 
cially when viewing angle constraints are available, can 
help narrow down the high altitude emission zone. A 
larger pulsar sample also allows a study of evolution of 
the 7-ray beaming and efficiency with pulsar age: the pul- 
sars seen with Fermi, not including millisecond pulsars, 
span ages from 10 3 to 2 x 10 6 yr (Abdo et al. 2009i), with 
hopes to extend that to larger r (lower E) as observa- 
tions continue. Analysis of the population of pulsars with 
interpulses (Weltevrede & Johnston 2008b) and radio po- 
larization data (Tauris & Manchester 1998) have given 
hints that the magnetic inclination is larger for young 
pulsars and decreases with age. Thus we may even probe 
evolution of magnetic alignment during pulsar spindown. 
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TABLE 1 

This table summarises the results of the timing and spectral analysis of PSR J1048— 5832 and PSR J2229+6114. Statistical 

AND SYSTEMATICS ERRORS ARE REPORTED. 



Analysis Parameters PSR J1048-5832 PSR J2229+6449 

Timing results Number of pulsed 7 rays 933 ± 93 1365 ± 97 

Peak position (<j>) 0.15 ± 0.01 ± 0.0001 (PI) 0.49 ± 0.01 ± 0.001 

0.57 ± 0.01 ± 0.0001 (P2) 

Peak separation (A) 0.42 ± 0.01 ± 0.0001 

Peak FWHM 0.06 ± 0.01 (PI) 0.23 ± 0.03 



0.10 ± 0.02 (P2) 



Spectral results a F (10~ 7 cm^s" 1 ) 


2.19 ± 0.22 ± 0.32 


3.77 ± 0.22 ± 0.44 


b F E (IQ- 11 erg cm^s- 1 ) 


19.4 ± 1.0 ± 3.1 


23.7 ± 0.7 ± 2.5 


r 


1.38 ± 0.06 ± 0.12 


1.85 ± 0.06 ± 0.10 


C E C (GeV) 


2. 3+04 ± 0.3 


3.6+°;^ ± 0.6 
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